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Abstract: We directly compare the concordance ACDM model to the inhomogeneous 
matter-only alternative represented by LTB void models. To achieve a "democratic" con- 
frontation we explore ALTB models with non- vanishing cosmological constant and perform 
a global likelihood analysis in the parameter space of cosmological constant and void radius. 
In our analysis we carefully consider SNe, Hubble constant, CMB and BAO measurements, 
marginalizing over spectral index, age of the universe and background curvature. We find 
that the ACDM model is not the only possibility compatible with the observations, and 
that a matter-only void model is a viable alternative to the concordance model only if the 
BAO constraints are relaxed. Moreover, we will show that the areas of the parameter space 
which give a good fit to the observations are always disconnected with the result that a 
small local void does not significantly affect the parameter extraction for ACDM models. 
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1. Introduction 

In the past 15 years it has been extensively studied (see for example [1, 2, 3, 4, 5, 6, 7, 8, 
9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]) how an observer inside a matter-only 
spherical void expanding faster than the background sees apparent acceleration. This effect 
is easy to understand: our cosmological observables are confined to the light cone and hence 
temporal changes can be associated with spatial changes along photon geodesies. In the 
present case, "faster expansion now than before" is simply replaced by "faster expansion 
here than there" . This is why a void model can mimic the effect of a cosmological constant 
if it extends to the point in space/time where the dark energy becomes subdominant. A 
typical scenario that can mimic the late-time acceleration of the concordance ACDM model 
consists of a deep void extending for 1-3 Gpc. 

Void models feature the heavy fine tuning of the observer's position which, to be 
consistent with the CMB dipole, has to be very close to the void center (few percents 



- 1 - 



[23, 24], se also [25]), thus leading to a violation of the Copernican principle. However, it 
is not clear how to regard this fine tuning as compared, for example, to the fine tuning 
of the cosmological constant or of other dark energy sources, and a pragmatic approach is 
just to rely on the observational data in order to find the most satisfying model. 

In the present work we adopt the latter pragmatic approach and to "democratically" 
find the best-fit model we study void models with cosmological constant. More specifically 
we will consider the so-called ALTB models, which are spherically symmetric dust Lemaitre- 
Tolman-Bondi (LTB) models with non- vanishing cosmological constant [26, 27, 28]. We 
will run likelihood analyses over two parameters: one (the density parameter Qa) pertinent 
to the ACDM model and the other (the void radius vq) pertinent to the LTB model. 1 Our 
parameter space will therefore contain both the ACDM model and the matter-only LTB 
model as delimiting axes, allowing a direct confrontation between the two opposite alter- 
natives and overcoming some of the ambiguities in comparing the two models. We remind 
that the LTB model is specified not by free parameters, but by free functions and therefore 
the number of degrees of freedom to be adopted in computing the reduced x 2 , which is 
usually used in model ranking, is necessarily a somewhat subjective choice. Our analysis 
will also show if the data favor a mixed scenario of a void together with a cosmological 
constant, which is a possibility not yet fully explored. 

It is worth pointing out that the study of void models represents just a small (even 
if admittedly the most successful) part of the research devoted to understand how large- 
scale inhomogeneities affect the observations, a phenomenon collectively referred to as 
backreaction (see [29]). Weak backreaction studies focus on the observational properties 
of the universe and include, besides the void models here examined, also swiss-cheese 
[30, 31, 32, 33, 34, 35, 36, 37], onion [38, 39] and meatball [40, 41, 42, 43] models. Strong 
backreaction studies, on the other hand, address the question of how the cosmological 
background reacts to the nonlinear structure formation, see for example [44, 45, 46, 47, 48, 
49, 50, 51, 52, 53, 54, 55, 56, 57]. 

This paper is organized as follows. In Section 2 we will introduce the specific ALTB 
model we will use, and in Section 3 we will discuss how to confront the ALTB model against 
the observations. In Section 4 we will show our results. We will see that the ACDM model 
is not the only possibility compatible with the observations, and that a matter-only void 
model is a viable alternative to the concordance model only if the BAO constraints are 
relaxed. Moreover, we will show that the areas of the parameter space which give a good 
fit to the observations are always disconnected. In particular, we will see that a local void 
does not affect the parameter extraction for ACDM models if its radius is smaller than 
1-2 Gpc. Finally, we will give our conclusions in Section 5 and we will discuss a possible 
degeneracy between 5q and ft,K,out in Appendix A. 

Together with this paper we release the Mathematica package LLTB 1.0, which is 
available at the address turbogl.org/LLTB.html. 



1 Also the void depth strongly characterizes a void model and we will keep it fixed to the value required 
by a good fit to the SNe data, see Section 4. 
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2. The ALTB model 



We will now quickly review the conventional LTB formalism. The general picture is of 
a local LTB void exactly matched to the chosen Friedmann-Lemaitre-Robertson- Walker 
(FLRW) background. For more details see, for example, Ref. [8]. 

2.1 Basic formalism 

The line element of the spherically symmetric LTB model can be written as (c = 1): 

ds 2 = -dt 2 + Y ' 2 ^\ dr 2 + Y 2 (r, t)(d6 2 + sin 2 Odcp 2 ) , (2.1) 
1 — k(r) 

where Y(r, t) is the scale function, the prime denotes derivation with respect to the coor- 
dinate radius r and the arbitrary function k(r) represents the local curvature. The FLRW 
solution is recovered by setting Y(r, t) — > a(t)r 2 and k(r) — >■ kr 2 throughout the equa- 
tions. Note that in the LTB space the transverse expansion rate Ht = Y/Y will generally 
differ from the longitudinal expansion rate Hl = Y'/Y'. For later use we also define the 
background expansion rate by H out (t) = a{t)/a{t). 

The dynamics of the model is governed by the following equation [28] : 

Y 2 _F(r) 8nG k(r) 

Y2 - ~yr + ~%~ PA ~y 2 ' ^ ' 

where the dot denotes derivation with respect to the coordinate time t and p\ = A/8nG is 
the energy density associated with the cosmological constant. The arbitrary function F(r) 
(actually a constant of integration) represents the effective gravitating mass and is related 
to the local dust energy density pM(r,t) through the equation F' = 8ttG pmY'Y 2 . It is 
useful to rewrite Eq. (2.2) in the following, more familiar form 



H 2 (r,t) = H 2 (r) 



O w ( ) I 3 + n A (r) + /s ; r - ( 



(2.3) 



Y(r,t) J v ' y ' \Y(r,t) 

where i?o( r ) = HT(r,to), Yq{t) = Y(r,to) and the (present-day) density parameters are 

w™m- (2 - 5) 

fii-(r) S 1 - S)„M - n A (r) = - H 2 ( ^ (r) • (2-6) 

We point out that these density parameters have their usual meaning (ratios of energy 
densities over the critical density) only within the homogeneous patches (origin included) 
where Hl = Ht, whereas generally the local critical density is given by p c = (H T + 
2HxHl) /8ttG. For later use we define also the present-day expansion rate at the observer's 
position -ffo.in = H$(r = 0) and the background present-day expansion rate i?o,out = 
Ho(r 3> rg), which clearly coincides with H ont (to). 
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Eq. (2.3) can be used to determine the age of the universe at a radial coordinate r: 

1 

t -t B (r) = —V- / , dx (2.7) 

We choose to constrain the models by requiring a simultaneous big bang, i.e., by setting 
t' B (r) = 0. Furthermore, we set the moment of the initial singularity at t = 0, so that to 
is the actual age of the universe. Simultaneous big bang excludes decaying modes which 
would be strongly in contradiction with the inflationary paradigm [58]. 

2.2 Specific model and its parameter space 

LTB models feature three arbitrary functions. Within the present formalism they are taken 
as ^m(0 (or F(r)), ie(r) (or k(r)) and Yo(r). One of these is but an expression of the 
gauge freedom, which we fix by setting Yo(r) = r and, consequently, a(to) = 1 where a(t) is 
the scale factor of the background FLRW model. Another free function, tsir), was already 
set by the simultaneous big bang condition and so we are left with specifying the matter 
profile Q.m(t). We will adopt the parameterization of Ref. [11] (constrained GBH) which 
reads: 

o o /o r> * 1 — tanh(r — ro/2Ar) 

fiM(r) = ^M,out + (fiAf,in " fiAf,out) ! + tanh(ro/2Ar) > ( 2 " 8 ) 

where the parameters ro and Ar characterize respectively size and steepness of the density 
profile, while &M,m and ^M,out are the matter density parameters at the observer's position 
and in the FLRW background outside the void. 

The precise form of the density profile should not be essential. The LTB void models 
(and void models in general) depend crucially on only two physical parameters: the void 
depth <5q = (fiM,in ~~ ^M,out)/^M,out which gives the jump AH in the expansion rate 
required to mimic the measured acceleration, and the void size ro which sets the redshift 
extension of the void. Any strong dependence on the precise shape of the profile (in our 
case the steepness Ar) would in fact signal fine tuning. 

We remind that our local LTB void is matched to the outside background FLRW 
model (where it is Y = a r) and so the local matter density profile exhibits a compensating 
overdense shell surrounding the central underdensity. We also point out that, as one can 
see from Eq. (2.4), Qm,ui is calculated with respect to the local critical density, i.e., using 
the higher local expansion rate. Therefore the void depth function 5n does not exactly 
correspond to the actual matter contrast 5m', for the parameter range we use in this work 
we typically find that 8m ~ <fo + 0.1. 

Summarizing, the background FLRW model will be specified by the parameters ^A.out, 
^■K,o\xt and to, while the void is modelled by 8q, ro and Ar. Moreover, we will leave the 
spectral index n s free and so the overall parameter space will be seven dimensional. 

3. Cosmological data analysis 

In this Section we will carefully explain how to compare the ALTB model predictions 
for an observer at its center with supernovae (SNe), Hubble constant, cosmic microwave 
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background (CMB) and baryon acoustic oscillations (BAO) observations. Before examining 
each of these datasets individually, we need to develop the formalism of the effective model 
necessary to confront CMB and BAO data. 

The LTB metric is matched to the background FLRW model at some radius r. Because 
the LTB solution corresponds to a pressureless dust source we enforce, similarly to Ref. [14], 
the matching to happen at a redshift at which radiation is still negligible. The value z = 100 
satisfies this requirement and allows the void to be large enough so as not to constrain 
artificially the parameter space. The relevant physics leading to CMB and BAO features 
occurs at redshifts much greater than z and so it is possible to describe the relative light 
cone by means of an effective FLRW model. One could be tempted to use the Einstein- 
de Sitter model for fitting because at z cosmological constant and curvature are locally 
negligible. The curvature, however, does not enter only the Friedmann equation, but also 
the metric with the global effect of changing, for example, the angular diameter distance. 
Because, differently from Ref. [14, 22], we will consider background models with nonzero 
curvature, we will use, similarly to Ref. [20], the FLRW background itself as the effective 
model. 

The calculation proceeds as follows. After finding the present-day expansion rate 
profile Ho(r) which satisfies Eq. (2.7) (this must be done numerically because no analytic 
solution can be found for i?o( r ) when A ^ 0), we solve the dynamics with Eq. (2.3). We 
then solve the light cone using the following geodesic equations: 



dt _ Y'(r,t) 

dz (l + z)Y'{r,t) ' ' } 

dr = y/1 - k(r) 

dz (l + z)Y'(r,t) ' K ' ' 

with initial conditions t(0) = to, r(0) = 0. As it is clear from the metric (2.1), the angular 
diameter and luminosity distance are then simply given by: 



d A (z) = Y(r(z),t(z)), (3.3) 
d L (z) = (1 + z) 2 d A (z) . (3.4) 

The effective model and the LTB model have to give the same angular diameter distance d A 
at the matching spacetime point (f, t) = (r(z),t(z)). To achieve this we place the effective 
observer at r = and we solve the light cone for the background FLRW model (ratio of 
Eqs. (3.1-3.2) in the FLRW limit): 

# = --^2=, (3-5) 
dr Vl - fcr 2 

with initial condition t(f) = i. The angular diameter distances now coincide because 
outside the void the scale function Y matches the FLRW scale factor : a(t) f = Y(f,t) = d A - 
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Finally, the effective metric is simply the background FLRW model at time t(r = 0) 
to, e ff- To be explicit, it is specified by the following parameters: 



Ho,eS = #out(*0,eff) , (3.6) 
l + Z , 
1 + Z e g 



l + Z 

To, e s = 1 - ?o , (3.7) 



tt 7 , cff = 2.469 • 10- 5 h~£ j , (3.8) 

n RteS = 1.692 7ieff , (3.9) 

^A,eff = "77^ ^A,out , (3.10) 

-"O.eff 

n _ ff o,out a 2 (to) , . 

^M,eff = 1 - ^ii,eff - ^A.cff - &>K,efi , (3-12) 

O s ,eff = 0.0226 h^i (^) • (3-13) 

The first two lines define the effective Hubble constant and the effective CMB tempera- 
ture. The CMB temperature measured by the real observer is set to To = 2.725 K and 
z c s = a (to,es)/ a (t) — 1- The dimensionless Hubble constant h e s is, as usual, defined by 
Ho,eft = 100/ieff km s _1 Mpc -1 . Eqs. (3.8-3.13) specify the effective density parameters; 
the modifications are due to the different time and CMB temperature experienced by the 
effective observer. Eq. (3.13) comes from requiring a fixed baryon-photon number density 
ratio [59]: the void observer's value is fixed to the WMAP7 [60] result of Q B h 2 = 0.0226, 
which is also compatible with Big Bang Nucleosynthesis constraints [61]. 

We will typically find that the effective CMB temperature is close to the observed one 
To, e s — Tq. This is mainly due to the fact that we are considering compensated voids, whose 
metric matches the outside background metric. The consequence is that for observables 
outside the LTB patch changes in redshift are small (see for example Ref. [31]). While this 
choice seems natural if one wants to recover a homogeneous universe on very large scales, 
it constrains the LTB model and its ability to fit the observables. This is similarly true 
for the simultaneous big bang condition we imposed with Eq. (2.7). We will develop these 
ideas in our Conclusions. 

3.1 Supernovae observations 

We will use the recent Union2 SNe Compilation [62] , which consists 557 type la supernovae 
in the redshift range z = 0.015—1.4. The predicted magnitudes are related to the luminosity 
distance of Eq. (3.4) by: 

m(z) = 51og 10 d L (z)/10pc, (3.14) 
and so the likelihood analysis is based on the x 2 function: 

^ = g H-"M + ri' , (315) 
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where the index i labels the Union2 entries. The parameter fi is an unknown offset sum of 
the SNe absolute magnitudes, of /c-corrections and other possible systematics. As usual, 
we marginalize the likelihood L' SNe = exp(— x'sNe/^) over Z 2 ' ^SNe = / dfJ-L' SNe , leading 
to a new marginalized x 2 function: 

XSNIa = S2-^, (3.16) 

where we dropped a cosmology-independent normalizing constant and the auxiliary S n is 
defined by: 



Note that since \jl is degenerate with log 10 -ffo,in we are effectively marginalizing also 
over the observer's Hubble constant. This is particularly important within the analysis 
of ALTB models where we can have models with very small but deep voids for which the 
value of the Hubble constant at the center is not as relevant in fitting the SNe as it is for 
large voids. 

3.2 Hubble constant 

The Hubble constant is obtained by measuring cosmological standard candles mostly within 
a distance of roughly 200 Mpc [63]. We therefore compute the local Hubble constant H\ oc 
by averaging the expansion rate profile within a sphere of radius r\ oc = 200 Mpc: 

/•Hoc / 

H Xoc = J flo(r) 4vrr 2 dr/ (4vr/3 rf oc ) , (3.18) 

where we used the fact that Yq(t) = r and we neglected the here irrelevant curvature. 

This averaging is meaningful in general, but it is of particular importance for an 
exploration of ALTB models. Indeed one can immagine a mixed scenario consisting of the 
concordance model plus a very small {tq <C 200 Mpc) but deep void. The unaveraged 
expansion rate at the very origin i^o,in would then differ significantly from the averaged 
H\ oc which better reflects the observable data being close to the background expansion rate 

#0,out- 

The present-day determination of H\ oc is converging but different groups still disagree 
on the precise value. We will consider the following two results from Ref. [64] and Ref. [65]: 

H S06 = 62.3 ± 5.2 km s _1 Mpc _1 Sandage et al. 2006 , (3.19) 

H R09 = 74.2 ± 3.6 km s^Mpc -1 Riess et al. 2009 . (3.20) 

We will perform the likelihood analysis separately for the two results, defining the respective 
X 2 functions as: 

2 _ (H SO q - H loc ) 2 

XS06 - 2 ' \ 6 - 11 ) 

°S06 

, 2 _ (#R09 ~ #loc) 2 00 x 

XR09 — 2 • (6.ZZ) 

°R09 
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3.3 Cosmic microwave background 

We will use positions and amplitudes of peaks and troughs in the CMB spectrum to con- 
strain the ALTB models. The location of peaks and troughs can be parametrized as [66]: 

l m = (m - <f> m ) I a , (3.23) 

where the peaks are labelled by integer values of m and the troughs by half-integer values, 
the quantity 4> m is a phase-shift parameter determined by pre-recombination physics and 
I A is the acoustic scale which is given by: 

l, -.r <'<* , >? + * , > . (3.24, 

r s 

where angular diameter distance, sound horizon r* and redshift z* at recombination are 
computed using the effective model of Eqs. (3.6-3.13). In our analysis we will consider 
the position of the first, second, third peak and of the first trough. We will compute the 
corresponding phases <j>i, 0i.5, 02 and 4>s using the accurate analytical fits of Ref. [67]. We 
will also consider the relative heights of second and third peak relative to the first one, H 2 
and H3, for which we can use the fits of Ref. [66]. 

We constrain the ALTB model by means of the following x 2 function built from the 
six quantities explained above: 



2 (Zi - Ziyw) 2 (h.5 - h.b,W7) 2 (h - h,W7 



,2 



XWMAP — I ^2 ! l~ a 2 

h h.5 h 

(h ~ k,W7) 2 (H 2 - H 2: w7) 2 (H 3 - H 3:Wr ) 2 
a? a% a% 

where the values marked with W7 correspond to the best-fit WMAP7 spectrum [60]. The 
error has to take into account both the experimental error and the error due to the fit. We 
thought it reasonable to use a a of 1% for the position of the first peak and of 3% for the 
other quantities. We checked the accuracy of the fits by comparing their predictions with 
the ones of CAMB [68] within our parameter space. 

We would like to stress that the quantity Xwmap ^ s a function of the spectral index n s . 
In particular, because it is computed through analytical fits, it is an analytical function 
over which we can easily integrate. We will see in Section 4 the importance of marginalizing 
over the spectral index. 



3.4 Baryon acoustic oscillations 

A baryon acoustic peak is detected from SDSS and 2dFGRS galaxy catalogues at redshifts 
0.2 and 0.35 [69]. The observed quantity is A# 2 Az, where A# is the angle that the acoustic 
scale subtends on the sky in the transverse direction at the observed redshift, and Az is 
the redhisft interval corresponding to the acoustic scale in the radial direction. In order 
to test void models against the BAO results, we have to understand how to compute the 
quantity A9 2 Az in a void model. This procedure is explained in detail in Ref. [20]. 
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First we need the comoving acoustic scale at the drag epoch rf ras , which is defined 
as the time at which the baryons are released from the Compton drag of the photons. 
Generally, in an inhomogeneous universe the drag epoch corresponds to different times in 
different locations: tj, = td{r). However, because we constrained the LTB model to have a 
simultaneous big bang (see Section 2.1), the universe is almost homogeneous at the drag 
epoch and so the r dependence in t^{r) is very weak. rf rag can then be computed using 
the effective metric and the fitting formulas of Ref. [70]. The latter, however, are tuned for 
models with a CMB temperature of To = 2.725 K and so we have to correct the relevant 
quantities with the appropriate powers of Tq/Tq^r. 

Because in the LTB metric transverse and longitudinal expansion rates differ (Ht ^ 
Hl), the acoustic scale is stretched differently in the transverse and radial directions. The 
correct scaling is given by the following formulas: 

l ^ = m$M>^- < 3 - 26 > 

where Zg rag (z) is the proper acoustic scale at the drag epoch for the shell r(z). As for 
td(r), the z dependence in ^ rag (z) is very weak because our early universe was close to 
homogeneity and we can simply write: 

drag 

e«=^, (3-28) 

where z<i is calculated using the effective metric. 

It is also interesting to see how the calculation proceeds in the general case where one 
cannot assume a weak r dependence in td(r). The final result for A6 2 Az is indeed formally 
the same. We start by noting (this can be seen from Eq. (2.3)), that radial shells evolve as 
independent FLRW models, and we assume that Az is small enough so that the spherical 
shells at r(z) and r(z + Az) can be described by the same FLRW model at the time t(z), 
where z is the redshift corresponding to the BAO observations. Next we evolve this FLRW 
model to find an observer that would see the event at time t(z) with the same redshift 
z and we use this observer's cosmological parameters to calculate rf rag and z& with the 
formulas in Ref. [70], and then use Eq. (3.28) to find Zf ag . 

The quantities AO and Az can then be written as: 

Az(z) = l L (z) (1 + z) H L {r{z),t{z)) , (3.30) 

where in the last line we have used Eq. (3.2) in the limit k(r), Az <C 1. Finally, we need 
the general relation: 

i I -n i ,( Y\r(z),t(z))Y 2 (r(z),t(z)) , n 
+ d ~ { + ' \Y'(r(z)Mr)) Y*(r(z)Mr)) ) ' ( ' 
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which shows that the volume element at r(z) has diluted by a factor (1 + £<z)/(l + z) 
between td(r) and t(z), i.e., that the density at drag epoch is correctly scaled. Lastly, we 
can construct the quantity (z~ 1 A6 2 Az) 1 / 3 = rf rag /<iy (z) by identifying 



where we have used Eq. (3.31). Note that the redshift Zd has canceled in the final expression; 
this would not be the case if one were to consider Az and AO separately. 

We will compare the void models with the results of Ref. [69]: rf ag /dy(0.2) = 0.1980± 
0.0058 and rs rag /dy(0.35) = 0.1094 ± 0.0033, by computing the \ 2 from the multivariate 
Gaussian likelihood. That is with 



4. Results 

As summarized at the end of Section 2, the parameter space of our ALTB model is seven 
dimensional and consists of the parameters fiA )0 utj ^K,out, to and n s for the background 
FLRW model and Sq, ro and Ar for the void. In this analysis we want to focus on 
(possibly) matter-only voids and so we set 5q = —0.9 (which corresponds to a matter 
contrast of 5m ~ —0.81) and Ar = 0.35 ro, because this choice always guarantees a good 
fit to the SNe. We will explore in Appendix A different values of <5q. In particular, we will 
show that there is a degeneracy between Sq and FlK,out as far as a good fit to the SNe is 
concerned. 

We are left with five parameters. We choose ^A,out and ro as the parameters with 
respect to which to plot the likelihood surfaces. These parameters indeed best characterize 
the two opposite alternatives of the ACDM and matter-only void models, which delimit 
the parameter space. We stress that this analysis will compare the two alternatives in an 
unbiased manner, thus overcoming some of the ambiguities in comparing the two models. 
We remind that the LTB model is specified not by free parameters, but by free functions 
and therefore the number of degrees of freedom to be adopted in computing the reduced \ 2 ■> 
which is usually used in model ranking, is necessarily a somewhat subjective choice. 

Only three parameters are now left to specify: n s , to and fi/f,out- We will plot in the 




1/3 



(3.32) 



X 2 bao = X T V- 1 X, 



(3.33) 



where 
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Figure 1: 1, 2 and 3a confidence level contours on tq and r?A,out for the likelihood of Eq. (4.1) 
for the ALTB model described in the text with n s — 0.96, to = 13.7 Gyr and ^K,out — 0. The 
four smaller panels on the left show the contours for the independent likelihoods per observable, 
while the larger panel on the right shows the contours for the combined likelihood. In the panel 
representing measurements of the local Hubble constant, the results relative to i?so6 °f Eq. (3.19) 
are shown as filled contours, while the ones relative to i?R09 of Eq. (3.20) are shown as lines. The 
same labelling holds for the panel relative to the combined observables. In the panel relative to 
CMB contraints we also show confidence levels for the likelihood of Eq. (4.2) marginalized over n s 
(dot-dashed contours). See Section 4.1 for a discussion. 
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(4.3) 


Fig 


. 4 


£(ro,^A,out) 


= J dn s dto dO,K, out L(ra, ^A.out, ^K,ont, to, n s ) , 


(4.4) 


Fig 


. 5 


£(?~0,^A,out, 


n s ) = j dto d9,R, out L(r , QA,out> ^Jf,out> to,n s ) , 


(4.5) 



where n s , to and ^x,out are some fixed values. Throughout the analysis we will adopt the 
prior ro < 3.5 Gpc: larger values of the void radius would result in a void whose light 
cone extends till times at which radiation cannot be neglected (see the discussion about 
the effective metric in Section 3). We will now show in Section 4.1 an illustrative example. 
We will then present our main results in Sections 4.2-4.5. 

4.1 Illustrative example 

We start by choosing n s = 0.96, to = 13-7 Gyr and ^K,out = 0, which are the WMAP7 
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values corresponding to the standard ACDM cosmology [60]. The results of the likelihood 
analysis are shown in Fig. 1. The four small panels on the left show the 1, 2 and 3a 
confidence level contours for the observables examined in the previous Section. We will 
now discuss them in detail. 

The SNe panel shows the known result that a pure-matter void model is as successful as 
the concordance ACDM model in fitting SNe observations. Moreover, the present analysis 
of ALTB models shows that a mixed scenario of a ACDM model together with a large deep 
void is not favoured. For smaller contrasts (as explained above we have fixed 5q = —0.9) the 
contours corresponding to the void model move towards the ones of the concordance model. 
However, the area of the parameter space with small vq and ^A.out (roughly delimited by 
the line joining concordance model and matter-only void) is excluded. For example, a 
shallow void with 5q = —0.4 and ^A,out = 0.6 still requires a radius of at least 1 Gpc. 

The "Local H" panel shows the constraints coming from local measurements of the 
Hubble constant. We remind the reader that we have fixed the age of the universe to and 
so the background expansion rate -£fo,out increases along the ^A,out axis. To understand 
the shape of the contours it is useful to look back at Eq. (3.18). If the void is large, then 
.f/ioc will be higher than the background value i^o.out because of the AH jump caused 
by Sq. If, however, the void radius is small, then the averaging of Eq. (3.18) gives back 
the lower -ffo.out value. The latter trend, together with the fact that i^o.out increases with 
^A,out ; shows that the -£fi oc relative to a point in the parameter space can also be obtained 
for smaller voids and higher ^A.out- For the parameters chosen, a matter-only model is 
consistent with the -ffso6 of Eq. (3.19), while it is not with the i/R.09 of Eq. (3.20). 

The CMB panel shows the constraints coming from the CMB spectrum. We show both 
the confidence levels relative to n s = 0.96 (filled contours and Eq. (4.1)) and the ones for the 
likelihood marginalized over n s G [0.8, 1.2] (dot-dashed contours and Eq. (4.2)). The fits we 
are using are valid in the latter range, which includes the values relevant for the likelihood 
analysis. For the parameters chosen the concordance ACDM model, possibly with a large 
local void, is favoured in both cases. In the non marginalized case the contours depend 
weakly on tq because the effect of the void on the angular diameter distance is generally 
small for observables outside the LTB patch. The dependence on r$ becomes stronger for 
very large voids which extend deep into the past light cone. 

The BAO panel shows the constraints coming from the measurements of the acoustic 
oscillations in the matter power spectrum. For this observable, the concordance model does 
not give a good fit, even though it is within the 3a confidence level contour. It is interesting 
to see that a la confidence level fit is achieved by a ACDM model with rU.out = 0.55 — 0.75 
together with a void of radius r$ = 0.8 — 2.7 Gpc. For the chosen to an d ^x,out) a pure- 
matter void model is excluded. We point out that by fine tuning the void profile beyond the 
Ansatz of Eq. (2.8) it is possible to fit the BAO data without cosmological constant [20]. 

Finally, the large panel on the right in Fig. 1 shows the combined likelihood (for 
n s = 0.96) with the result that the standard concordance model is favoured against a void 
model of large radius of any type for these parameters. As shown in the inset, a small local 
void of few tens Mpc is not, however, excluded. Note in particular that the confidence 
level contours are perpendicular to the ^A,out axis and so a small local void does not bias 
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Figure 2: 1, 2 and 3er confidence level contours on ro and f2A,out for the likelihood of Eq. (4.2) for 
the ALTB model described in the text marginalized over n s and with to = 13.7 Gyr and f2.R-,out = 0. 
Labelling as in Fig. 1. The right panel shows the combined likelihood for all observables, while the 
BAO constraints are excluded in the left panel. See Section 4.2 for a discussion. 



the parameter extraction within the framework of ACDM models. See Ref. [71] for an 
alternative analysis. 

It is important to point out that the previous results depend crucially on the values 
chosen for n s , to and ^K,out- We decide, therefore, to marginalize the total likelihood over 
the latter parameters in order to analyze the data with as few as possible priors. We show 
these results in the next Sections, while the findings of this Section should be considered 
as illustrative ones. We stress that in this way we do not need either the reduced \ 2 or the 
minimum \ 2 to compare the ACDM with void models. 

4.2 Marginalizing over n s 

In this Section we will consider models with to = 13.7 Gyr and ST2A',out = but of any 
spectral index, that is, we marginalize the total likelihood over n s as shown in Eq. (4.2). 
Numerically we consider a range n s 6 [0.8,1.2], which includes the values relevant for 
the likelihood analysis, whose results we plot in Fig. 2. The right panel corresponds to 
a combined likelihood for all observables with the result that the concordance model is 
favoured at 3<7 confidence level against a void model of large radius of any type. 

In deriving the BAO contraints we assumed that the BAO scale is comoving. Whether 
this really is a good approximation is still unclear [72]. However, it is difficult to improve 
the treatment of the BAO feature because the perturbation theory in the LTB background 
is not yet fully understood. Therefore, in the left panel of Fig. 2 we repeated the analysis 
excluding the BAO contraints. In the case of the likelihood with i?so6 of Eq. (3.19), the 
la confidence level now contains, besides the concordance model, an (almost) pure-matter 
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Figure 3: 1, 2 and 3er confidence level contours on ro and f2A,out for the likelihood of Eq. (4.3) for 
the ALTB model described in the text marginalized over n s and to and with f^A'.out = 0. Labelling 
and definitions as in Fig. 2. See Section 4.3 for a discussion. 



void model. In the case of the likelihood with -fzR09 of Eq. (3.20) only the 3a confidence 
level contains the pure-matter void model. 

It is clear from our results that the areas of the parameter space which give a good fit 
to observations are always disconnected. Moreover, as pointed out in the previous Section, 
the confidence level contours relative to the concordance model are perpendicular to the 
^A,out axis. These two facts combined together show that a local void does not affect the 
parameter extraction for the ACDM model if its radius is smaller than 1-2 Gpc. The same 
conclusions hold also for the results of the next Sections. 

4.3 Marginalizing over n s and to 

In this Section we will consider models with ClK,out = but of any age and spectral 
index, that is, we marginalize the total likelihood over to and n s as shown in Eq. (4.3). 
Observations of stellar globular clusters are reported to result in a lower bound of to = H-2 
Gyr [73]. This bound, however, uses a ACDM template in order to relate the age of the 
oldest globular cluster in the Milky Way (10.4 Gyr) to the age of the universe. Relaxing the 
assumption of the background model we obtain a slighter lower bound of to — H Gyr. We 
then take to = 16.7 Gyr as upper bound, which we have found to be numerically equivalent 
to an unbounded upper limit as the total likelihood is already essentially zero. 

The results are shown in Fig. 3 where, similarly to the previous Section, we include 
(right panel) and exclude (left panel) BAO constraints in the analysis. With the BAO 
constraints the concordance model is favoured at 1 and 2a confidence level against a void 
model of large radius of any type. At 3a confidence level, however, a mixed scenario 
with a void of radius ro = 3.25 — 3.5 Gpc and f2A,out = 0.25 — 0.45 is not excluded. 
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Figure 4: 1,2 and 3er confidence level contours on r and f2A,out for the likelihood of Eq. (4.4) 
for the ALTB model described in the text marginalized over n s , to and f2_R- i0 ut- Labelling and 
definitions as in Fig. 2. See Section 4.4 for a discussion. 



Without the BAO constraints the results are similar to the ones of the previous Section if 
Hsqq is used. The contours are instead different in the case of the likelihood with .ffRog, 
which at la confidence level now has a mixed scenario with a void of large radius and 
^A.out = 0.15 — 0.30. The concordance model is always within the la confidence level 
contour. 

4.4 Marginalizing over Tisi t$ and ^Kout 

In this Section we will consider models of any spectral index, age and background curvature, 
that is, we marginalize the total likelihood over n s , to and ^K,ont as shown in Eq. (4.4). 
As in the previous Sections, n s £ [0.8,1.2] and to 6 [H> 16.7] Gyr. The background 
curvature range we consider is 0_R- j0 ut £ [—0.25, 0.16] which includes the values relevant 
for the likelihood analysis. The results are shown in Fig. 4 where again we include (right 
panel) and exclude (left panel) BAO constraints in the analysis. 

Without the BAO contraints, the la confidence level fully includes the matter-only 
void model if Hso6 is used, while the concordance model is always included. If -ffRog is 
used, the pure-matter void model is included only by the 3a confidence level. With BAO 
constraints included the la confidence level always contains the concordance model but 
not the matter-only void model. However, the 3a (filled contours relative to Hsoq) or 2-3a 
(empty contours relative to i^Rog) confidence levels include a mixed scenario of a large void 
with f2A,out ~ 0.1 — 0.5. 

We would like to point out that without the BAO constraints the matter-only void 
model can give as good a minimum \ 2 as the concordance model if we use i?so6 for the 
local Hubble parameter. For example, a pure-matter void model with to = 15.45 Gyr and 
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Figure 5: 1, 2 and 3er confidence level contours on r and f2A,out for the likelihood of Eq. (4.5) 
for the ALTB model described in the text marginalized over to and ^lx,out and with n s — 0.96. 
Labelling and definitions as in Fig. 2. See Section 4.5 for a discussion. 



^if,out = —0.1 gives a minimum \ 2 of 559.8, while the concordance model a \ 2 °f 559.4, 
where in both cases we marginalized over n s . 

We would like to comment that, within our modelling, a curved background can have 
two different interpretations. fix,out can indeed be relative to the global universe or it 
could just describe a patch bigger than the LTB void. In the latter case one could imagine 
the presence of many patches of different £Ik, out s ° that on average the universe is still flat. 

4.5 Marginalizing over to and f^out 

Finally, in this Section we will consider models with n s = 0.96 but of any age and back- 
ground curvature, that is, we marginalize the total likelihood over to and ^Kout as shown 
in Eq. (4.5). As in the previous Sections, to £ [H> 16.7] Gyr and &K,out £ [—0.25, 0.16]. 
The results are shown in Fig. 5 where again we include (right panel) and exclude (left 
panel) BAO constraints in the analysis. 

With the BAO constraints the concordance model is favoured at 3a confidence level 
against a void model of large radius of any type. The same is true without the BAO 
constraints if ^ro9 is used. In the case of -Hso6) instead, pure- matter void models are 
within the 2 and 3a confidence levels, with the la confidence level favouring models of 
large radius and f^A.out ~ — 0.2. 

5. Conclusions 

In the present paper we attempted to "democratically" confront the concordance ACDM 
model with the inhomogeneous matter-only alternative represented by the LTB void mod- 
els. To this end we performed a likelihood analysis of ALTB models in a seven-dimensional 



- 16 - 



parameter space describing a local void embedded in a possibly curved FLRW background. 
We plotted likelihood surfaces with respect to the background cosmological constant f2A iOU t 
and void radius tq. These two parameters best characterize indeed the limiting cases of 
ACDM and pure-matter void models, which are recovered as delimiting axes of the pa- 
rameter space. We set the void depth to the value required by a good fit to the SNe 
observations. 

In order to analyze the data with as few priors as possible, we marginalized the total 
likelihood over the spectral index n s , the age of the universe to an d the background cur- 
vature f^out- See Eqs. (4.1-4.5) for a quick summary and legend. In this way, it is not 
necessary to compute the reduced x 2 , or even its minimum, thus overcoming some of the 
ambiguities in comparing the two models. We remind that the LTB model is specified not 
by free parameters, but by free functions and therefore the choice of the number of degrees 
of freedom to be adopted in computing the reduced x 2 > which is usually used in model 
ranking, is necessarily made on a somewhat subjective basis. 

We confronted the ALTB model with SNe, Hubble constant, CMB and BAO observa- 
tions and we found that the concordance model is not the only possibility compatible with 
the observations. In particular, by allowing a nonzero curvature, we found that a large void 
with r^A.out ~ 0.2 — 0.4 lies within the 2a confidence level contour, while a pure-matter 
model is excluded. 

However, these results depend on the precise treatment of the BAO constraints. If we 
do not include BAO contraints in the analysis - we remind that perturbation theory in an 
LTB background is not thoroughly studied yet - a matter-only model can be as successful 
as the concordance model. We stress that our modelling of the void depends crucially only 
on the void radius and depth, which are the two main physical quantities describing an 
underdensity, while by a more specified tuning of the void profile it is possible to fit the 
BAO data without cosmological constant [20]. 

There is a number of ways one could improve the analysis of the present paper and 
obtain better results for the void model. First, one could scan a larger parameter space. 
A possibility regarding the void parameters could be to let free the density contrast (see 
Appendix A) and one regarding the background parameters could be to consider a running 
spectral index. Second, one could extend the modeling of the void; we found interesting 
four directions. 

1. One could drop the simultaneous big bang condition we imposed with Eq. (2.7). 
Simultaneous big bang excludes decaying modes which would be strongly in contra- 
diction with the inflationary paradigm [58], but Gpc scale voids are anyway at odds 
with the standard scenario. Without simultaneous big bang one would have a new 
free function which would likely improve the agreement with the observations. See 
for example Ref. [11]. 

2. By tuning the density profile one can find better agreement with the BAO obser- 
vations and with the local measurements of the Hubble rate. Void models already 
feature the heavy tuning of the observer position, which has to be close to the center 
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in order not to give a too large dipole. An extra fine tuning on the density profile 
could be treated in the same pragmatic way. See for example Ref. [20]. 

3. Uncompensated voids can give sizeable redshift effects for observables outside the 
LTB patch, while this is not generally true for compensated voids. In particular it 
would be interesting to study very large voids with vq > 3.5 Gpc. In order to perform 
such an analysis, however, it would be necessary to generalize the ALTB formalism 
in order to include the radiation: very large voids extend deeply into the past light 
cone where radiation cannot be neglected. See for example Ref. [14] . 

4. The possibility that the radiation density is not uniform is also interesting. As in 
the previous point, it requires a proper modeling of the radiation content. See for 
example Ref. [21]. 

By exploiting these extensions it may be possible to successfully fit all the present-day 
cosmological observables. In particular, the recent claim [74] that void models give a 
too large kinetic Sunyaev-Zel'dovich effect could be accommodated within a more general 
modeling. 

Finally, we have seen that the areas of the parameter space which give a good fit to the 
observations are always disconnected. Our analysis has indeed shown that a local void does 
not affect the parameter extraction for ACDM models if its radius is smaller than 1-2 Gpc. 

Together with this paper we release the Mathematica package LLTB 1.0, which is 
available at the address turbogl.org/LLTB.html. We also would like to recommend an 
interested reader to visit the same webpage, where pre-compiled Mathematica notebooks 
with animations for exploring a wide range of the to and &K,out parameters are available. 
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A. Degeneracy between 5q and fi/^out 

In the analysis of Section 4 we set 8q = —0.9 (which corresponds to a matter contrast of 
5m ~ —0.81) because it always guarantees a good fit to the SNe. If FIk, out / 0, however, 
the latter is not a necessary requirement as far as a good fit to the SNe is concerned. The 
reason is that the background curvature can account for part of the necessary jump in AH 
demanded to mimic the observed acceleration. 

In order to explore this possibility we have plotted in Fig. 6 two examples of models 
with a lower contrast that give good fits to the SNe. The parameter values are the same of 
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Figure 6: 1, 2 and 3cr confidence level contours on ro and r^out for the likelihood of Eq. (4.1). For 
the panels on the left it is f2_R- i0 ut = 0.2 and <5ji — —0.85 (which corresponds to a matter contrast 
of w —0.75). For the panels on the right it is £lfc,out = 0.4 and Sq — —0.8 (which corresponds 
to a matter contrast of « —0.7). In both cases it is n s — 0.96 and to = 13-7 Gyr. Labelling and 
definitions as for the four left panels in Fig. 1. See Appendix A for a discussion. 



Fig 1 with the difference that, instead of Ox,out = and 5n = —0.9, it is £lK,out = 0.2 and 
<5r> = —0.85 (which corresponds to a matter contrast of ~ —0.75) for the panels on the left 
and r^out = 0.4 and 5q = —0.8 (which corresponds to a matter contrast of rj —0.7) for the 
panels on the right. As one can see, a lower contrast needs a higher background curvature, 
thus confirming a degeneracy between Sq and ^x.out- Moreover, an open background 
mimics an underdensity and so also the void radius is smaller. 

We have seen that an open background allows to have a good SNe fit with a shallower 
and smaller void, surely a desirable scenario. To conclude on the viability of this setup 
we have to look at the other observables. About local H constraints, by comparing Fig. 1 
and Fig. 6 we see that the confidence level contours are basically unaffected as far as pure- 
matter models are concerned. The reason is that higher rix,out gives a higher i?o,out which 
is (al least partially) compensated by the smaller AH of the shallower void. About the 
remaining observables, while the confidence level contours relative to the BAO observations 
change in favor of shallower voids, a good fit to the CMB needs a close background with 
the result that the models of Fig. 6 are strongly ruled out and the same conclusion holds 
if we consider different values for the age of the universe. 

Concluding, relaxing the prior we have used for the matter contrast could slightly 
widen the contours towards smaller values of the void radius, but not strongly affect the 
results of Section 4. Finally, we would like to stress that the particular contrast needed for 
a good fit to the SNe depends on the particular density profile chosen. 
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